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Abstract 

Jet production is studied in the Breit frame in deep-inelastic positron-proton scattering over 
a large range of four-momentum transfers 5 < < 15 000 GeV^ and transverse jet ener- 
gies 7 < Et < 60 GeV. The analysis is based on data corresponding to an integrated lumi- 
nosity of £int — 33 pb~^ taken in the years 1995-1997 with the HI detector at HERA at a 
center-of-mass energy ^/s = 300 GeV. Dijet and inclusive jet cross sections are measured 
multi-differentially using k± and angular ordered jet algorithms. The results are compared 
to the predictions of perturbative QCD calculations in next-to-leading order in the strong 
coupling constant a^. QCD fits are performed in which and the gluon density in the 
proton are determined separately. The gluon density is found to be in good agreement with 
results obtained in other analyses using data from different processes. The strong coupling 
constant is determined to be as(M^) = 0.1186 ib 0.0059. In addition an analysis of the 
data in which both and the gluon density are determined simultaneously is presented. 
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1 Introduction 



Deep-inelastic lepton-proton scattering (DIS) experiments have played a fundamental role in 
establishing Quantum-Chromodynamics (QCD) as the theory of the strong interaction and in 
the understanding of the structure of the proton. The lepton inclusive DIS cross section is 
directly sensitive to the quark densities in the proton, but gives only indirect information on 
the gluon content and on the strong coupling constant via scaling violations of the structure 
functions. The production rates of events in which the final state contains more than one hard 
jet (besides the proton remnant) are, however, observables which are directly sensitive to both 
as and the gluon density in the proton. These multi-jet cross sections can thus be used to test the 
predictions of perturbative QCD (pQCD) and allow a direct determination of as and the gluon 
density [jl]]. 

The large center-of-mass energy ^/s of 300 GeV at HERA allows multi-jet production in 
DIS to be studied over large regions of phase space. In this paper we present comprehensive 
measurements of jet production in the range of four-momentum transfers squared 5 < < 
15 000 GeV^. Using four different jet algorithms we study multi-differential distributions of the 
dijet and the inclusive jet cross sections to which the predictions of pQCD in next-to-leading 
order in as are compared. We identify those observables for which theoretical predictions have 
small uncertainties and perform QCD analyses of the jet data in which we determine the value 
of a^. A consistent determination of the gluon density in the proton, together with the quark 
densities is obtained in a simultaneous fit which additionally includes data on the inclusive DIS 
cross section. In the last step an analysis of the data in which both as and the gluon density are 
determined simultaneously is presented. 

This paper is organized as follows. In section 2 we give a short description of the theoret- 
ical framework and motivate the choice of the jet variables to be measured. The experimental 
environment and details of the measurement procedure are described in section 3 and the multi- 
differential jet cross sections are presented in section 4. Finally, in section 5 we introduce the 
theoretical assumptions and the methods which are used in the QCD analysis and present the 
results of the QCD fits. Numerical values of the results are given as tables in the appendix. 



2 Jet Production in Deep-Inelastic Scattering 
2.1 Jet Variables and the Breit Frame 

The inclusive neutral current cross section in deep-inelastic lepton-proton scattering is described 
in lowest order perturbation theory as the scattering of the lepton off a quark in the proton via 
the exchange of a virtual gauge boson (7, (according to Fig. [1] (a)). The kinematics of the 
reaction are given by the four-momentum transfer squared Q^, the Bjorken scaling variable xb] 
and the inelasticity y defined as 



{I - I'f 



p-q 



(1) 



2p ■ q 



y = 



p ■ I 
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where / (/') and p are the four-momenta of the initial (final) state lepton and proton, respectively. 
When particle masses are neglected the kinematic variables are related to the lepton-proton 
center-of-mass energy ^/s by s xbj y = ■ The variable xbj is in the leading order approxima- 
tion identical with the longitudinal momentum fraction x of the proton which is carried by the 
parton specified by the parton density functions, hereafter referred to as the struck partonQ. 

Multi-jet production in DIS is described by the QCD-Compton and the boson-gluon fusion 
processes. Due to the latter contribution multi-jet cross sections are directly sensitive to the 
gluon density in the proton. Examples of leading order diagrams of both processes are shown 
in Fig. [1] (b) and (c). Both diagrams contribute to the cross section which depends explicitly on 
as- Variables that characterize features of the multi-jet final state are the invariant mass Mjj, the 
partonic scaling variable Xy and the variable ^ defined as 

e = a:Bj(l + -^) and x,, = ^ . (2) 

In the leading order picture (when the final state partons are identified with jets) the invariant 
dijet mass Afjj is equal to the center-of-mass energy of the boson-parton reaction. In this ap- 
proximation the fractional momentum x of the struck parton is given by the variable ^ which 
becomes much larger than the Bjorken scaling variable x^^ if Mjj is large. The partonic scaling 
variable x^ specifies the fractional momentum of the incoming parton seen by the boson. 

Studies of the dynamics of multi-jet production are preferably performed in the Breit frame 
where the virtual boson interacts head-on with the proton [||]. The Breit frame is defined by 
2xBjP + g = 0, where p and q are the momenta of the proton and the exchanged boson, respec- 
tively. The positive 2;-axis is chosen to be the proton direction. In the lowest order process, at 
C(a°), the quark from the proton is back- scattered into the negative ^-direction and no trans- 
verse energy is produced^. The appearance of jets with large transverse energies Et can only 
be explained by hard QCD processes whose contribution is at least of 0[oi^ relative to the 

'in this paper the Bjorken scahng variable xbj is always written with a subscript to distinguish it from the 
proton momentum fraction x which appears in the formulae of the proton's parton densities. While the former is 
an observable quantity, the latter is only defined in a theoretical framework within a given factorization scheme. 

^By "transverse" we refer to the component perpendicular to the z-axis. The transverse energy is defined as 
Et = EsmB. The polar angle Q is defined with respect to the proton direction in both the laboratory frame and 
the Breit frame. Throughout the paper "transverse energy" always refers to transverse energies in the Breit frame. 
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gluon-photon center-of-mass frame 
X p = - q i ^ / 




Breit frame 



XgjP = -1/2q 




Breit,2 



z-boost 



Pt 



2^V 



2ti' 



Figure 2: A boson-gluon fusion event in deep-inelastic scattering in the boson-gluon center- 
of-mass frame (left) and in the Breit frame (right). The frames are related to each other by a 
longitudinal boost along the z-direction. 

inclusive DIS cross section. The hardness of the QCD process is specified by Et which is the 
physical scale at which e.g. hard gluon radiation is resolved. 

In the leading order approximation the Breit frame is related to the boson-parton rest frame 
by a longitudinal boost along the 2;-direction (see Fig. |^). The polar scattering angle of the 
jets in the boson-parton center-of-mass frame is directly related to the pseudorapidity[| i]' of 
the jets. In leading order approximation the value of r]' is equal to half the difference of the 
jet pseudorapidities r/Breit in the Breit frame rj' = llr^Breit.i — ^Breit,2|- Since the transverse jet 
energy Et is invariant under longitudinal boosts along the 2-axis Et is identical in both frames. 

2.2 Jet Definitions 

The comparison of the properties of high multiplicity hadronic final states observed in the ex- 
periment to those in perturbative calculations involving only a small number of partons requires 
the definition of infrared- and coUinear-safe jet observables. While the properties of different jet 
observables depend on the exact definition of the jets the physical interpretation of experimen- 
tal results must, however, not depend on details of the jet definition if theory is to be claimed 
successful. 

In this analysis we use four jet clustering algorithms which successively recombine particles 
into jets. All jet algorithms are applied in the Breit frame to the final state particlesQ excluding 
the scattered lepton. They can be grouped into two pairs of inclusive and exclusive jet algo- 
rithms, each pair consisting of one kj_ ordered and one angular ordered algorithm. In the kj_ 
(angular) ordered algorithms pairs of particles are clustered in the order of increasing relative 
transverse momenta k± (increasing angles) between the particles. The exclusive jet definitions 
assign each particle explicitly to a hard jet or to the proton remnant, while for the inclusive 

^The pseudorapidity is defined as 77 = — ln(tan 6*72) where 9 is the polar angle. Positive values of ij correspond 
to particle momenta pointing into the proton hemisphere. For massless particles, differences in pseudorapidity are 
invariant under longitudinal boosts. 

""^"Particle" refers in this paper either to an energy deposit or a track in the detector, to a parton in a perturbative 
calculation or to a hadron (i.e. any particle produced in the hadronization process including soft photons and leptons 
from secondary decays). All particles are treated as massless by a redefinition of the energy (E = \p\). 
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jet definitions not all particles are necessarily assigned to hard jets. The following four jet 
algorithms are used: 

• the exclusive k± ordered algorithm as proposed in [|^]. 

• the exclusive angular ordered algorithm (Cambridge algorithm) as proposed in [Q] and 
modified for DIS to consider the proton remnant as a particle of infinite momentum along 
the positive z-axis, following the approach used in [|3|]. The exact definition is taken 
from [|5|]. 

• the inclusive k_L ordered algorithm as proposed in [7]] . 

• the inclusive angular ordered algorithm (Aachen algorithm) as proposed in In 
analogy to the changes from the exclusive k± algorithm to the Cambridge algorithm, the 
inclusive algorithm has been modified to obtain an inclusive jet algorithm with angular 
ordering. 

The recombination of particles in the exclusive jet algorithms is made in the .E-scheme 
(addition of four-vectors) resulting in massive jets. To maintain invariance under longitudinal 
boosts for the inclusive jet definitions the Et recombination scheme [Q] is used in which the 
resulting jets are massless. 

In the exclusive jet definitions the clustering procedure is stopped when the distances i/ij = 
^"^ defined between all pairs of jets and between all jets and the proton remnant are above 
some value ?/cut, where k'^-j = 2mm{Ef, Ej){l — cos%) and 5 is a reference scale. In our 
analysis we set S"^ = 100 GeV^ and ?/cut = 1 to have the final jets separated by k^ij > 10 GeV. 
The inclusive jet algorithms are independent of an explicit stopping criterion in the cluster- 
ing procedure and hard jet selection cuts have to be applied afterwards. These algorithms are 
defined by a radius parameter Rq which we set to i?o = 1 as suggested in 

2.3 Phase Space 

The kinematic region in which the analysis is performed is defined by the kinematic variables y 
and g2 

0.2 <y < 0.6 and 5 < < 15 000 GeV^ . (3) 

The lower limit on y has been chosen to exclude the kinematic region of large sbj where jets 
are predominantly produced in the forward direction, i.e. at the edge of the detector acceptance. 
The upper limit on y ensures large energies of the scattered lepton. 

The jet finding is performed using the jet algorithms introduced above. We restrict the 
jet phase space to the angular range in which jets can be well measured in the HI detector. 
Therefore the four-vectors of the jets defined in the Breit frame are boosted to the laboratory 
frame where we then apply the pseudorapidity cut 

-1 < %t,lab < 2.5 . (4) 
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In this kinematic range double-differential jet cross sections are measured as a function of Q^, 
Et (inclusive jet cross section) and Et = |(-E'ri + Et2) (dijet cross section) using the various 
jet algorithms mentioned above. In addition the dependences of the dijet cross section on the 
dijet variables Mjj, Xp and i]' as introduced in section ^J] as well as on the pseudorapidity 
of the forward jet r/forw, lab in the laboratory frame are measured. In all cases inclusive dijet 
cross sections, i.e. cross sections to produce two or more jets within the angular acceptance are 
measured. The jet variables are calculated from the two jets with highest transverse energy (the 
jets are labeled in the order of descending Et). 

In the measurement of the dijet cross section care has to be taken to avoid regions at the 
boundary of phase space which are sensitive to soft gluon emissions where perturbative calcu- 
lations in fixed order are not able to make reliable predictions. The exclusive jet algorithms 
avoid these regions due to the cut on the variable k±. For the inclusive jet definitions additional 
selection cuts have to be chosen appropriately. The dijet cross section defined by a symmetric 
cut on the transverse energy of the jets -Et, 1,2 > -f'T.min is infrared sensitive [[TOp. This prob- 
lem can be avoided by an additional, substantially harder cut on, for example either a) the sum 
Et,i + Et,2, b) Et,i or c) the invariant dijet mass Afjj. When cuts are chosen to obtain cross 
sections of similar size in all of a), b) and c) above, the next-to-leading order corrections are 
largest in b) and hadronization corrections are largest in c). The smallest next-to-leading order 
corrections and hadronization effects are seen for scenario a). For the inclusive jet algorithms 
we therefore require 

Et, 1,2 > 5 GeV and Et,i + Et,2 > 17 GeV . (5) 



2.4 QCD Predictions of Jet Cross Sections 

While leading order (LO) calculations can predict the order of magnitude and the rough features 
of an observable, reliable quantitative predictions require the perturbative calculations to be 
performed (at least) to next-to-leading order (NLO). The NLO calculations of the jet cross 
sections used in this analysis are performed in the MS scheme for five massless quark flavors 
using the program DISENT0 [ |TT1 ] which has been tested in [ [T2| ] and found to agree with the 



program DISASTER-i-i- [|T3]] in the kinematic region of interest. 

Perturbative fixed order calculations beyond leading order can give reliable quantitative pre- 
dictions for observables with small sensitivity to multiple emission effects and non-perturbative 
contributions. They fail, however, to predict details of the structure of multi-particle final states 
as observed in the experiment. A complementary approach to describe these properties of the 
hadronic final state is used in parton cascade models. Starting from the leading order ma- 
trix elements, subsequent emissions are calculated based on soft and coUinear approximations. 
There exist two different approaches in which parton emissions are either described by a parton 
shower model (HERWIG 0, LEPTO [0 and RAPGAP [[T|]) or by a dipole cascade (ARI- 



ADNE [JTTp). These parton cascade models can be matched to phenomenological models of the 



hadronization process. The HERWIG event generator uses the cluster fragmentation model QlSl ] 



while in LEPTO, RAPGAP and ARIADNE the Lund string model M is implemented. The 



programs HERWIG, LEPTO, RAPGAP and ARIADNE are used in the present measurement 
^We have modified DISENT to include the running of the electromagnetic coupling constant. 
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Figure 3: The predictions of (a) the hadronization corrections to the dijet cross section for 
different jet definitions as a function of as obtained by HERWIG and (b) the next-to-leading 
order corrections to the dijet cross section as a function of for the inclusive k±_ algorithm 
using two different renormalization scales 



to provide event samples which are used in the correction procedure for the data. In the QCD 
analysis they are used to estimate the size of the hadronization corrections to the perturbative 
jet cross sections. 

Higher order QED corrections can change the size of the cross section and also modify the 
event topology. Especially hard photon radiation may strongly influence the reconstruction of 
the event kinematics and thereby the boost vector to the Breit frame. Corrections from real 
photon emissions from the lepton and virtual corrections at the leptonic vertex are included in 
the program HERACLES [jg^] which is directly interfaced to RAPGAR An interface to LEPTO 
and ARIADNE is provided by the program DJANGO (grp. 

Safe predictions can only be expected for observables for which perturbative higher-order 
corrections and non-perturbative (hadronization) corrections are small. Detailed investigations 
on properties of the NLO cross sections and the size and the uncertainties of the hadronization 
corrections to the observables under study have been performed in [Bl 1^, P2p. The hadronization 



corrections predicted from the different models are in good agreement and have small sensitivity 
to model parameters. The hadronization corrections (5hadr. are displayed in Fig. |^ (a) as a func- 
tion of for aU jet algorithms used. They are defined as 5hadr. = (cThadron - o"parton)/o-parton 
where ciparton (o"hadron) IS the jct cross section before (after) hadronization. The hadronization 
corrections are generally smaller for the inclusive jet algorithms than for the exclusive ones and 
smaller for the ordered algorithms when compared to those with angular ordering. Hence 
the inclusive kx_ algorithm shows the smallest corrections, acceptable even down to very low 
Q"^ values. Having the smallest hadronization corrections, the inclusive k±_ algorithm is thus 
the best choice for a jet definition. The other jet algorithms will, however, still be used to 
demonstrate the consistency of the results. 

An indication of the possible size of perturbative higher-order contributions is given by the 
size of the NLO corrections or the renormalization and factorization scale dependence of an 
observable. For the inclusive ki_ algorithm the NLO corrections to the dijet cross section are 
displayed in Fig. ^ (b). Shown is the A;-factor, defined as the ratio of the NLO and the LO 



6 



predictions, for two different choices of the renormalization scale (/i^ = Et, Q). Towards low 
Q'^ the NLO corrections become large, especially for the choice fir = Q. Reasonably small k- 
factors (k < 1.4) are only seen at > 150 GeV^ where Q'^ and are of similar size such that 
terms oc \n{E^/Q'^) are small. The renormalization scale dependence is seen to be correlated 
with the NLO correction i.e. large at small Q"^. The factorization scale dependence is below 2% 
over the whole phase space (not shown). These studies suggest that a QCD analysis of jet cross 
sections, involving the determination of and the gluon density, should be performed at large 
values ofQ'^. 



3 Experimental Technique 

The analysis is based on data taken in the years 1995-1997 with the HI detector at HERA 
in which positrons with energies of E^ = 27.5 GeV collided with protons with energies of 
Ep = 820 GeV. 



3.1 HI Detector 

A detailed description of the HI detector can be found elsewhere . Here we briefly introduce 



the detector components most relevant for this analysis. 

In the polar angular range 4° < 6' < 154° the electromagnetic and hadronic energy is 



measured by the Liquid Argon (LAr) calorimeter [ |24| ] with full azimuthal coverage. The LAr 
calorimeter consists of an electromagnetic section (20—30 radiation lengths) with lead absorbers 
and a hadronic section with steel absorbers. The total depth of both sections varies between 4.5 
and 8 interaction lengths. Test beam measurements of the LAr calorimeter modules have shown 



an energy resolution of cte/E ^ 0.12/ ^yE [GeV] © 0.01 for electrons [|2|] and as/E 



0.50/ a/-E [ GeV] © 0.02 for charged pions after software weighting 

In the backward direction (153° < 9 < 177°) energy is detected by a lead-fiber calorime- 



ter, SR^CAL [27]. It consists of an electromagnetic section with a depth of 28 radiation 



lengths in which the scattered positron is measured with an energy resolution of ^e/E 



0.071/a/-E' [GeV] © 0.010. It is complemented by a hadronic section to yield a total depth 
of two interaction length. 

Charged particle tracks are measured in two concentric jet drift chamber modules (CJC), 
covering the polar angular range 25° < 9 < 165°. A forward tracking detector covers 7° < 
9 < 25° and consists of drift chambers with alternating planes of parallel wires and others with 
wires in the radial direction. A backward drift chamber BDC improves the identification of the 
scattered positron in the SR^CAL calorimeter. The calorimeters and the tracking chambers are 
surrounded by a superconducting solenoid providing a uniform magnetic field of 1.15 T parallel 
to the beam axis in the tracking region. 

The luminosity is measured using the Bethe-Heitler process ep e^p. The final state 
positron and photon are detected in calorimeters situated close to the beam pipe at distances of 
33 m and 103 m from the interaction point in the positron beam direction. 
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3.2 Event Selection 



Neutral current DIS events are triggered and identified by the detection of the scattered positron 
as a compact electromagnetic cluster. The data set is divided into two subsamples in which 
the positron is detected either in the SPACAL (5 < < 70 GeV^) or in the LAr calorimeter 
(150 < < 15 000 GeV^) with uniform acceptance over the range 0.2 < y < 0.6. These re- 
gions are labeled "low Q^" and "high Q^" throughout the text. The low Q'^ and high samples 
correspond to integrated luminosities of Ci^t — 21 pb^^ and £int — 33 pb^^, respectively^. 

At low the positron is reconstructed as the highest electromagnetic energy cluster in the 
SPACAL, requiring an energy of > 10 GeV and a polar angle of 156° < 9^ < 176°. The 
positron selection at high closely follows the procedure used in the recent measurement of 
the inclusive DIS cross section [pS]], requiring an electromagnetic cluster of E'^ > 12 GeV with 
a polar angle 9q < 153°. For 6'e > 35° the positron candidate is validated only if it can be 
associated with a reconstructed track, which points to the positron cluster. Fiducial cuts are 
applied to avoid the boundary regions between the calorimeter modules in the z and (f) (i.e. 
azimuthal) directions. The events in the low and high samples are triggered by demanding a 
localized energy deposition together with loose track requirements. The trigger efficiencies for 
the final jet event samples are above 98%. 

In both samples the reconstructed z-coordinate of the event vertex is required to be within 
±35 cm of its nominal position. The hadronic final state is reconstructed from a combination of 
low momentum tracks (pr < 2 GeV) in the central jet chamber and energy deposits measured in 
the LAr calorimeter and in the SPACAL according to the prescription in From momentum 
conservation the sum J2i^^Pz) over all hadronic final state particles and the scattered positron 
is expected to be 2Ee = 55 GeV. This value is lowered in events in which particles escape 
undetected in the beam pipe in negative ^-direction. Photoproduction background and events 
with hard photon radiation coUinear to the positron beam are therefore suppressed by a cut on 
45 <J2iE-Pz)< 65 GeV. 

The event kinematics is determined from a redundant set of variables using the scattered 
positron and the hadronic final state. Using all hadronic final state particles h the variable 
^ = Ylihi^h ~ P^a) derived. The kinematic variables xbj, y and Q'^ are then reconstructed 
by the Electron-Sigma method [ ]29| ] 

^2 AJ7 JP' 2 -E'e^ sin^ 9e 2Ee 

<5eS = 4^e^e COS — X^E = 7Z T = ^ , ^ Z/E 

2 s|/e(1-|/s) S + E^(1 - cos6'e) 



L + _£/g(l — COSt^e) 

The jet algorithms are applied to the hadronic final state particles which are boosted to the Breit 
frame. The boost vector is determined from the variables |/es, Qes the azimuthal angle of 
the scattered positron. The transverse jet energy Et (or Et), the dijet mass Mjj and the variable 
7]' are calculated from the four- vectors of the jets. The variables ^ and Xp are reconstructed as 



6ec = a^eS H and Xp^rec = ^ with ?/h = TTFT • (7) 

yh S 4rec 2Ee 



*The low sample uses only data from the years 1996-1997. 
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These relations exploit partial cancellations in the hadronic energy measurement in Mj? and y^. 

The fraction of dijet events in the inclusive neutral current DIS event sample varies strongly 
with Q^, namely between ~ 1% (at = 5 GeV^) and ~ 20% (at = 5000 GeV^). Using 
the inclusive kj_ algorithm we have selected 1 1400 dijet events at low and 2855 dijet events 
at high Q^. The inclusive jet sample (measured only at high Q^) contains 10432 jets with 
Et > 7 GeV, from 7263 events. The size of photoproduction background has been estimated 
using two samples of photoproduction events generated by PYTHIA [ ]30| ] and PHOJET pl|]. 
The contribution to the distributions of the finally selected events is found to be negligible (i.e. 
below 1%) in all variables under study. 



3.3 Correction Procedure 

The data are corrected for effects of limited detector resolution and acceptance, as well as for 
inefficiencies of the selection and higher order QED corrections. The latter are dominated by 
real photon emissions from the positron (initial- and final-state radiation) and virtual corrections 
at the leptonic vertex, as included in the program HERACLES. No further corrections for effects 
due to the running of the electromagnetic coupling constant or non-perturbative processes (i.e. 
hadronization) are applied. 

To determine the correction functions the generators LEPTO, RAPGAP and ARIADNE (all 
interfaced to HERACLES) are used. For each generator two event samples are generated. The 
first sample, which includes QED corrections, is subjected to a detailed simulation of the HI 



detector based on GEANT [p2|]. The second event sample is generated under the same physics 
assumptions, but without QED corrections and without detector simulation. The correction 
functions are determined bin-wise for each observable as the ratio of its value in the second 
sample and its value in the first sample. This method can be used if migrations between different 
bins are small and properties of the simulated events are similar to those of the data. The 
absolute normalization of the generated cross sections is, however, arbitrary, since this cancels 
in the ratio. 

To test their applicability for the correction procedure, detailed comparisons have been made 
of the simulated event samples and the data for a multitude of jet distributions [Q]. None of the 
models can describe the magnitude of the jet cross section; especially at low large deviations 
are seen. However, all models give a reasonable description of the properties of the hadronic 
final state and of the properties of single jets and the dijet system including angular jet distribu- 
tions. In a previous publication [ p3| ] it has been shown that these event generators give a good 
description of the internal structure of jets. 

Based on the event simulation the bin sizes of the observables are chosen to match the 
resolution. The final bin purities and efficiencies are typically above 50% and migrations are 
sufficiently small to have small correlations between adjacent bins. The correction functions 
as determined by different event generators are in good agreement with each other and the 
absolute values deviate typically by less than 20% from unity. The final correction functions 
applied to the data are taken to be the average values from the different models. The difference 
between the average and the single values are quoted as the uncertainty induced by the model 
dependence which is subdivided in equal fractions into correlated and uncorrelated uncertainty 
between data points. This uncertainty is typically below 4%. 
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3.4 Experimental Uncertainties 



In addition to the model dependence of the correction function and the statistical uncertainties 
of the data and of the correction function several other sources of systematic experimental un- 
certainties are studied. They are given in the following, together with the typical change of 
the cross sections and a remark whether a particular uncertainty is treated as correlated or un- 
corrected between different data points. The latter classification closely follows the one used 
in pi]. 



The measurement of the integrated luminosity introduces an overall normalization uncer- 
tainty of ±1.5%; correlated. 

The hadronic energy scale of the LAr calorimeter is varied by ±4%; ±2% of the effect is 
considered to be correlated; typical change of the cross sections ±7.5%. 

The hadronic energy scale of the SPACAL is varied by ±7%; 
typical change of the cross sections < ±1%; uncorrected. 

The track momenta of the hadronic final state are varied by ±3%; 
typical change of the cross sections ±2.5%; uncorrected. 

The calibration of the positron energy in the SPACAL is varied by ±1%; 
typical change of the cross sections < ±2%; correlated. 



The positron calibration of the LAr calorimeter is treated as in [[280; a variation between 
±0.7% and ±3% is made, depending on the 2; -position of the energy cluster in the detec- 
tor, from which ±0.5% is considered to be correlated between different data points; the 
rest is treated as uncorrected; typical change of the cross sections ±4%. 

The positron polar angle is varied by ±2 mrad (±3 mrad) for positrons in the SPACAL 
(LAr calorimeter); typical change of the cross sections < ±2%; correlated. 

The positron azimuthal angle is varied by ±3 mrad; 
typical change of the cross sections < ±1%; uncorrected. 



The largest experimental uncertainty comes from the uncertainty of the energy scale of the LAr 
calorimeter. Since the uncertainties from all other sources are fairly small their determination is 
often subject to fluctuations. We therefore give a conservative estimate, by quoting the maximal 
(up- and downward) variation as the symmetric uncertainty. 

The statistical and the uncorrected systematic uncertainties are added in quadrature to ob- 
tain the total uncorrected uncertainty. The correlated contributions are kept separately and can 
thus be considered in a statistical analysis. To obtain the total uncertainty for each single data 
point, all contributions are added in quadrature. 
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Figure 4: The inclusive jet cross section as a function of the transverse jet energy in different 
regions of Q'^ for the inclusive k± algorithm (left) and for the Aachen algorithm (right). The 
data are compared to the perturbative QCD prediction in NLO with (dashed line) and without 
(solid line) hadronization corrections included. 

4 Experimental Results 

The measured cross sections, corrected for detector effects and effects of higher order QED, 
are presented as single- or double-differential distributions where the inner (outer) error bars 
represent the statistical (total) uncertainty of the data points. The results (defined in the phase 



space specified in section p3| ) are directly compared to the perturbative QCD predictions in 
NLO. The hadronization corrections 5hadr. have been estimated using the models described in 



section for which the predictions are in good agreement with each other. In all the figures 
shown the theoretical prediction "NLO (8) (l+5hadr.)" is derived from the NLO calculations with 
hadronization corrections determined using HERWIG. All NLO calculations are performed us- 
ing the parton density parameterizations CTEQ5M1 and a value of as{Mz) = 0.118. 
The renormalization scale is set to the transverse jet energy fXr = Et or in case of the dijet 
cross section to the average transverse energy Et- For the factorization scale a fixed valueQ of 
fif = -\/200 GeV, corresponding to the average Et of the jet sample, is chosen. 



^This slightly unusual procedure is motivated in section 5.2. A variation of /x/ in the range 6 < /z/ < 30 GeV 
changes the NLO results by less than 2%. 
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Figure 5: The ratio of the measured inclusive jet cross section and the theoretical prediction 
for the inclusive k± algorithm (top) and the Aachen algorithm (bottom). The uncertainty of the 
theoretical prediction is indicated by the band (the contributions from the renormalization and 
factorization scale dependence and the hadronization corrections are added in quadrature). 



4.1 Inclusive Jet Cross Section 



The inclusive jet cross section is measured at high in the Breit frame for both inclusive 
jet algorithms. The results are presented in Fig. ^ double-differentially as a function of the 
transverse jet energy in the Breit frame Et in different regions of Q^. The data for the inclusive 
algorithm (left) and for the Aachen algorithm (right) cover a range of transverse jet energies 
squared (49 < < 2500 GeV^) which is similar to the range of the four-momentum transfers 
squared (150 < < 5000 GeV^) of the event sample. The cross sections are of the same 
size for both jet algorithms and show a slightly harder Et spectrum with increasing Q^. The 
hadronization corrections are seen to be below 10% for both algorithms. The ratio of data 
and theoretical prediction is shown in Fig. |5|. Over the whole range of Et and the NLO 
calculation, corrected for hadronization effects, gives a good description of the data. 



4.2 Dijet Cross Section 

The inclusive dijet cross section is measured over the large range of four-momentum transfers 
squared 5 < < 15 000GeV^ using the inclusive jet algorithm. At high Q"^ additional 
measurements have been made using the three other jet algorithms. 

The dijet cross section for the inclusive k± algorithm is shown in Fig. ^ as a function of 
in the range 5 < < 15 000GeV^ for the central analysis cut (Eti + Et2 > 17 GeV) 
and an additional harder cut (Eti + Et2 > 40 GeV) on the sum of the transverse energies of 
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Figure 6: The dijet cross section measured with the inclusive k± algorithm as a function of 
for different cuts on the sum of the transverse jet energies. The data are compared to the pertur- 
bative QCD prediction in NLO (solid line), in LO (dotted line) and to a theoretical prediction 
where hadronization corrections are added to the NLO prediction (dashed line). 



the two jets with highest Ex- The dataQ are compared to the NLO prediction without and with 
hadronization corrections applied, as well as to a LO calculation. 

The hadronization corrections are small and increase only slightly towards lower Q^. The 
NLO prediction, including hadronization effects, gives a good description of the data over the 
large phase space in Et and and nicely models the reduced dependence observed for 
the higher Et data. Deviations at Q'^ ~ 10 000 GeV^ can be attributed to the neglect of Z° 
exchange in the calculation. 

The dijet cross section at high Q"^, measured using all four jet algorithms mentioned above, 
is shown in Fig. 0. A different dependence is observed for the inclusive and the exclusive 
algorithms which is a reflection of the different jet selection criteria and which is well repro- 
duced by the theory. While hadronization corrections have only a small effect for the inclusive 
jet algorithms, they lower the NLO predictions for the exclusive algorithms by up to 30% at 

= 150 GeV^. However, when these non-perturbative corrections are included the dijet cross 
sections are in all four cases well described by the theoretical curves. 

In the following more details of the dijet distributions are given. For these studies we restrict 
the phase space to < 5 000 GeV^ in order to avoid the region where contributions from Z° 

**A11 cross sections have been measured as bin-averaged cross sections and all but two are presented this way, 
the only exception being the presentation of the dependence in Figs. ^jand^. Here, to compare the data to the 
differential NLO prediction the points are presented at the bin-center, as determined using the NLO calculation. 
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Figure 7: The dijet cross section as a function of for four jet algorithms. The data are 
compared to the perturbative QCD prediction in NLO with (dashed line) and without (solid 
line) hadronization corrections included. 



exchange are sizable. We present the double-differential dijet cross section as a function of 



the variables Q 



Mjj, Et, 



^ and ?7forw, lab in Figs. || - [T3|. As for the distribution in 



Fig. 0, the results are compared to the perturbative QCD prediction in NLO with and without 
hadronization corrections included. In Fig. [TT] the contribution from gluon-induced processes 
is shown in addition. 

The distributions of the invariant dijet mass Mjj and the average transverse jet energy 
are shown in Fig. |S| covering a range of 15 < Mjj < 95 GeV and 8.5 < Et < 60 GeV. 
In both distributions we observe a harder spectrum towards larger Q^. The NLO prediction, 
including hadronization corrections, gives a good overall description, except at lowest Q'^ where 
it describes the shape, but not the magnitude, of the cross section. The increasing hardness of 
the Et distribution at higher Q'^ is also seen for the other jet algorithms in Fig. 1^. 



The distribution of the pseudorapidity rj' (as defined in section [23| ) is shown in Fig. [T^ for 
different regions of Et for the low and the high data. In both data sets the fraction of jets 
produced centrally in the dijet center-of-mass frame is observed to be larger at higher Et. 

The partonic scaling variable Xp is defined as the ratio of the Bjorken scaling variable xbj 



and the reconstructed parton momentum fraction ^. In the distribution shown in Fig. 11 (left) 



a strong variation of the Xp range is seen. Towards lower Q values ^ differs by up to three 
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Figure 8: The dijet cross section for the inclusive k± algorithm as a function of Mjj (left) and 
Et (right) in different regions of Q^. The data are compared to the perturbative QCD prediction 
in NLO with (dashed line) and without (solid line) hadronization corrections included. 



orders of magnitude from XBy At leading order the variable ^ represents the fraction of the 
proton momentum carried by the struck parton. The dijet cross section in bins of ^ is therefore 
directly proportional to the size of the parton densities at the parton momentum fraction x = ^. 
Fig. |Tl] (right) shows the ^ distribution in different regions of Q^. The dijet data are seen to be 
sensitive to partons with momentum fractions 0.004 < ^ < 0.3 which only increase slightly 
with increasing Q^. The ^ distribution is of special importance in the QCD analysis for the 
determination of the gluon density in the proton. Therefore we display the contribution from 
gluon induced processes to this distribution which varies strongly from ~ 80% at low to 
~ 40% at the highest Q^. Both, the ^ and the Xp distributions are well described by the NLO 
calculation over the whole range of Q^, independent of the fractional gluon contribution. 
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Figure 9: The dijet cross section as a function of the average transverse jet energy in the Breit 
frame in different regions of Q"^ for the inclusive k± algorithm (top left), the Aachen algo- 
rithm (top right), the exclusive k± algorithm (bottom left) and the Cambridge algorithm (bottom 
right). The data are compared to the perturbative QCD prediction in NLO with (dashed line) 
and without (solid line) hadronization corrections included. 
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Figure 10: The dijet cross section for the inclusive algorithm as a function of the variable r]' 
in different regions of Et at low Q'^ (left) and at high Q'^ (right). The data are compared to the 
perturbative QCD prediction in NLO with (dashed line) and without (solid line) hadronization 
corrections included. 



The ^ distribution is also presented for the other jet algorithms (Fig. [T^. While the distri- 
butions for the inclusive jet algorithms (incl. k± and Aachen) are already described by the NLO 
calculation (without hadronization corrections), large deviations are seen for the exclusive al- 
gorithms (excl. k± and Cambridge), especially at small ^ corresponding to small dijet masses. 
However, in this region hadronization corrections are very large for the exclusive algorithms. 
Within the estimated size of these corrections theory and data are consistent, except in those 
regions where the corrections are especially large. 

Fig. |13] finally shows the distribution of the forward jet r^forw.iab in the laboratory frame in 
different regions of Q^. While at larger the distribution is seen to decrease towards the cut 
value at r^forw.iab = 2.5, it is flatter at low Q^. The theoretical calculation gives a reasonable 
description of this angular distribution. In addition also the LO prediction is included. Although 
the NLO corrections become large in the forward region (i.e. at large ?7for„,iab) towards lower 
Q"^, the NLO calculation does describe the data remarkably well. Only at lowest Q"^ the NLO 
calculation clearly fails to describe the data, which is in agreement with the observations made 
in an earlier analysis of forward jet production [^. 



17 



CM 
> 

^ 2000 

Q. 

^ 
^ 500 

T3 

Q. 

X 
T3 



lo 100 

CO 

T3 



20 




2 




1 




0.2 



0.01 



5<Q < lOGeV' 



• HI data 
NLO CTEQ5M1 
NLO (gluon only) 



10<Q^<20 GeV^ 



20 < < 35 GeV^ 



35 < < 70 GeV^ pZf 



-150<Q^<200 GeV^ 



a. 



200<Q^<300 GeV^ 



300 < < 600 GeV^ 



600 < < 5000 GeV^ 



incl. kj^ algo 



J] 



10 



10 



-2 



10 



> 

^ 1000 

Q. 



o 

T3 

■o 

■D 



500 



200 


100 




20 




10 




2 




0.1 



5<Q <10GeV'' 



• HI data 
NLO GTEQ5M1 
NLO (gluon only) 



1 < < 20 GeV^ 



• J 



20 < < 35 GeV^ 



35 < < 70 GeV^ 



1 50 < < 200 GeV^ 




-200<O^<300 GeV^ 




300 < 0^ < 600 GeV^ 




600 < O^ < 5000 GeV^ 



incl. kj^ algo 



10-^ 10-^ 

k -XBj(1 +M,5/q2) 



Figure 11: The dijet cross section for the inclusive k± algorithm as a function of the variables 
Xp (left) and ^ (right). The perturbative QCD prediction in NLO (solid line) is compared to 
the measured dijet cross section. In addition the contribution from gluon induced processes is 
shown (dashed line). 



The perturbative NLO prediction gives a good description of the data for those observables 
for which NLO corrections and non-perturbative contributions are small. This agreement is 
seen in all regions of phase space, independent of whether they are dominated by the QCD- 
Compton or the boson-gluon fusion processes. For observables with not too large hadronization 
corrections the differences between the perturbative calculation and the data can be explained 
by the predictions of phenomenological hadronization models. In the kinematic region of 10 < 

< 70GeV^ theory still gives a good description of the data although NLO corrections 
become large. The theoretical calculations only fail at < 10 GcV^ where NLO corrections 
are largest (with /c-factors above two), such that contributions beyond NLO are expected to be 
sizable. 
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Figure 12: The dijet cross section as a function of the reconstructed parton momentum fraction 
^. The data are measured in different regions of Q'^ for the inclusive k± algorithm (top left), 
the Aachen algorithm (top right), the exclusive k± algorithm (bottom left) and the Cambridge 
algorithm (bottom right). The data are compared to the perturbative QCD prediction in NLO 
with (dashed line) and without (solid line) hadronization corrections included. 
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Figure 13: The dijet cross section for the inclusive k± algorithm as a function of the pseudo- 
rapidity of the forward jet in the laboratory frame. The data are compared to the perturbative 
QCD prediction in NLO (solid line), in LO (dotted line) and to a theoretical prediction where 
hadronization corrections are added to the NLO calculation (dashed line). 
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5 QCD Analysis 



The QCD predictions depend primarily on as and on the gluon and the quark density functions 
of the proton. In this section we present QCD analyses of the data in which we determine 
these parameters of the theory. We briefly discuss how different processes in DIS are directly 
sensitive to the different parameters and introduce the physical and technical assumptions with 
which the QCD fits are performed. 



5.1 Strategy 

In perturbative QCD (pQCD) the cross section of any process in deep-inelastic lepton-proton 
scattering can be written as a convolution of (process specific) perturbative coefficients Ca,n with 
(universal) parton density functions fa/p of the proton 



a 



y2 / dx Ca,„ (— ,/ir-,/i/) /a/p(a;,/i/) • (8) 

„ „ Jo 



The sum runs over all contributing parton flavors a (quarks and gluon) and all orders n con- 
sidered in the perturbative expansion. The integration is carried out over all fractional parton 
momenta x. The coefficients Ca,n are predicted by pQCD. They are currently known to next- 
to-leading order in the strong coupling constant for the inclusive DIS cross section (n = 0, 1) 



and for the dijet and the inclusive jet cross section (n = 1,2) [pq]. In the regions of sufficiently 
large transverse jet energies and not too large values of (Q^ < 5 000 GeV^) the effects of Z° 
exchange and of quark masses (for five quark flavors) can be neglected as shown in |Q] using 
the program MEPJET p7| ] . In this approximation the perturbative coefficients of the quarks for 
the inclusive DIS cross section and for the jet cross sections fulfill the relations 

Cu = Cc = Cu = Cg and ca = Cg = Cf, = = Cg = ci (9) 

in each order of a^. Therefore only three coefficients are independent and the cross section in 
(H) can be described by three independent parton density functions^ xG{x), xA(x) and xT.(x) 
with coefficients cq, ca and cs, which have to be defined such that 

Cgg{x) + ^Caiqaix) + qaix)) = cgG{x) + csS(x) + caA(x), (10) 

a 

where the sums run over all quark flavors a. These three parton density functions are chosen to 
be 

Gluon: xG{x) = xg(x), 

Sigma: xS(x) = x {qajx) + qajx)) , 

a 

Delta: x A{x) = x ^ el{qa{x) + qa{x)) , (11) 



'We do not explicitly display the dependence on the factorization scale. 



21 





LO 


NLO 


0"mcl. DIS 


xA{x) 
xG{x), xA{x) 


xA{x), xG{x) 
xG{x), xA{x), xT,{x) 



Table 1 : Overview of the parton density functions contributing to different cross sections in LO 
and NLO. 

where Ca denotes the electric charge of the corresponding quark. The three corresponding coef- 
ficients are given by linear combinations of the single flavor coefficients 

Cg = Cgluon Cs = ^ (4 Q - C„) Ca = 3 (c„ - Cd) ■ (12) 

At orders and 0(al) the contributions from different quark flavors are proportional 

to their electric charge squared (i.e. Cu = ^ Cd). Therefore the coefficient ce in ( [T^ vanishes and 
the only quark contributions to the cross sections come from xA(x). The gluon gives contribu- 
tions at order 0{al) and higher. starts to contribute at order 0{al) and does therefore 
not enter the inclusive DIS cross section to next-to-leading order. Table [1] gives an overview of 
the orders in which the parton densities contribute to the different processes (to NLO). xS(x) 
enters only the jet cross sections via the NLO corrections. At large the contributions from 
xS(x) are, however, small (4.5% for the dijet cross section at 150 < < 200 GeV^, decreas- 
ing to 2% at 600 < Q2 < 5000 GeV^). In the following the parameterization CTEQ5M1 is 
used to determine this contribution which is not regarded as a degree of freedom in the analysis. 
This is, however, only a weak assumption which will (due to the smallness of the contribution) 
not bias the result. 

With this approximation the inclusive DIS cross section and the jet cross section now depend 
on three quantities which will be determined in this analysis: a^, the gluon density xG{x) and 
the quark density xA(x). To demonstrate the basic sensitivity the leading order cross sections 
are written in the symbolic form 

inclusive DIS cross section: cTinci. dis A 

jet cross sections in DIS: a^ct oc ag-iccG + ca A) . (13) 

These relations make clear that in DIS a direct determination of either as or the gluon density 
can never be performed without considering the correlation with the other quantity. Three 
strategies are used in the QCD fits which differ by the amount of external information included 
in the analysis. 

1. Determination of from jet cross sections: using the jet cross sections measured one 
can determine assuming external knowledge on the parton distributions as provided 
by global data analyses. 

2. Consistent determination of the gluon density xG{x) and the quark density xA{x): in- 
cluding data on the inclusive DIS cross section, which are directly sensitive to the quarks 
only, and assuming the world average value of as{Mz) the information provided by the 
jet data can be used for a direct determination of the gluon density together with the quark 
densities via a simultaneous fit. 
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3. Simultaneous determination of a^, the gluon density xG{x) and the quark density xA(x): 
if the jet cross sections are measured in different phase space regions (ajet, cjet) with 
different sensitivity to the quark and the gluon contributions (i.e. where c'^/ c'^ ^ cg/c^) a 
simultaneous direct determination of all free parameters is possible when again additional 
inclusive DIS data are included. 



5.2 Fitting Technique 

A determination of theoretical parameters can only be performed in phase space regions where 
theoretical predictions are reliable. Although the perturbative NLO calculation gives a good 
description of the jet data down to = 10 GeV^, the QCD analysis is restricted to the region 
where NLO corrections are small (with A;-factors below 1.4), i.e. to the region of high (150 < 
< 5 000 GeV^). For the main analysis the jet cross sections measured for the inclusive k± 
algorithm are used for which hadronization corrections are smallest. Jet cross sections from 
other jet algorithms are used to test the stability of the results. The uncertainties of the jet data 
and their correlations are treated as described in section . 



In the second and in the third step of the analysis data on the (reduced) inclusive DIS cross 
section are included to exploit their sensitivity to the quark densities in the proton. A subsample 
is taken of the recently published measurement in the range 150 < < 1 000 GeV^. 
Since the present analysis uses the same experimental techniques the effects of the point to 
point correlated experimental uncertainties can be fully taken into account. 

The fit of the theoretical parameters is performed in a minimization using the program 
MINUIT . The definition of [ |39| ] fully takes into account all correlations of experimental 



and theoretical uncertainties. This definition has also been used in recent global data analy- 
ses pO| , |4T| ] and in a previous H 1 publication . The quoted uncertainties of the fit parameters 



are defined by the change of the parameter for which the x"^ of the fit is increased by one. 

In the fitting procedure the perturbative QCD predictions in NLO for the inclusive DIS cross 
section are directly compared to the data, while the NLO predictions for the jet cross sections 
are corrected for hadronization effects before they are compared to the jet data: 



HI ^ ^ ^NLO 

"^incl. DIS ^ ^ '^incl. DIS 



^hadron _ parton 



C^iet ^ ^ ^mt ■(l + 'Jhadr.) With dhadr " 



'jet ^jet V-^ I "hadr.; "vini i/jiadr. — parton 

'^jet 



The hadronization corrections are determined as described in section using the average 



value from the model predictions by HERWIG, LEPTO and ARIADNE. The uncertainty from 
the model and the parameter dependence of these predictions is always below 3% |^]. The 
uncertainty in the matching of the parton level (parton cascade and NLO calculation) is taken 
into account by increasing the quoted uncertainty in those regions where the hadronization 
corrections are large. In detail, the uncertainty of the hadronization correction for each bin 
of the jet cross sections is taken to be half the size of the correction, but at least 3%. This 
uncertainty is assumed to be correlated between the theoretical predictions for all data points. 

The renormalization scale fir in the NLO calculation is identified with the process specific 
hard scales in both processes. The inclusive DIS cross section is evaluated at /i^ = ^/O^ and 
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the jet cross sections are evaluated at /i^ = Ex (inclusive jet cross section) and /i^ = Et (dijet 
cross section). For the jet cross section an example is given of how the results change for an 
alternative choice, = a/Q^- "^^^ strong coupling constant as{fir) is parameterized in terms 
of its value at the scale /i^ = Mz using the numerical solution of the renormalization group 
equation in 4-loop accuracy]^ P3|, 



In principle the arguments invoked in the choice of the renormalization scale fir also apply 
to the factorization scale jif for the inclusive DIS cross section and for the jet cross section. 
However, a different choice is made for the following reasons. The different parton flavors 
have been combined into three independent parton density functions xG{x, fif), xA{x, fif) and 



/i/) (section pTTl) . These three parton densities are, however, only independent as long 
as no evolution between different scales /i/ is performed. The evolution of the gluon density is 
coupled to the evolution of /i/). Furthermore, since xA{x, fif) is not an eigenstate of the 
DGLAP evolution operators, the evolution requires its decomposition into a non-singlet and a 
singlet (i.e. a;S(x, /i/)). This introduces an additional dependence between the quark densities. 
To avoid mixing between the different parton densities the parton distributions are not evolved 
to different scales. Instead the perturbative calculations are carried out at a fixed value of the 
factorization scale jij = fiQ. The jet cross sections are sensitive to the parton distributions in 
the x-range 0.008 < s < 0.3 (see Fig. [11]). In this x-range the factorization scale dependence 
of the parton density functions is not large. In a next-to-leading order calculation the remaining 
fif dependence given by the DGLAP evolution equations is largely compensated by a corre- 
sponding term in the perturbative coefficients. The perturbative cross sections therefore depend 
only weakly on the choice of the factorization scale. The difference between using a fixed fac- 
torization scale fiQ and performing the full DGLAP evolution at a scale /x/ is of higher order 
in as than those considered. If the scale /i/ is close to the fixed scale these higher order terms 
which are proportional to h\(fif/ hq) are small. Therefore a fixed value of the factorization scale 
of the order of the average transverse jet energies in the dijet and the inclusive jet cross section 
fif = fiQ = V200 GeV ~ {Et) is used. The subsample of the (reduced) inclusive DIS cross 
section 150 < < 1 000 GeV^ has been chosen such that the four- momentum transfer is also 
of the same order of magnitude a/Q^ ^ jjQ = y/200 GeV. 

Both the renormalization and the factorization scale dependences of the cross sections each 
are considered as correlated theoretical uncertainties. Both scales are (separately) varied by a 
factor Xu around their nominal values un in the range x„ = 2| and the ratios ""'^loCj^^m-w) 
taken as the corresponding uncertainties. Together with the uncertainty from the hadronization 
corrections they constitute the quoted theoretical uncertainty of the fit results. 

During the minimization procedure in the fit the NLO calculations of the jet cross sec- 
tions have to be performed iteratively for different values of as{Mz) and for different parton 
density functions (the number of calculations used to obtain the present results and to study their 
stability is in the order of one million). Since standard computations of NLO jet cross sections 
are time consuming the method is used of pre-convoluting the perturbative coefficients with 
suitably defined functions which can then be folded with the parton densities and for a fast 
computation of the NLO cross section. 

"'it has been checked that in the range of scales considered in this analysis, 7 GcV < jir < Mz, the differences 
between the 2-, 3- and 4-loop solutions are always below 3 per mille. 
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5.3 Determination of 



As a first step the QCD predictions are fitted to the jet cross sections using parameterizations for 
the parton distributions from global fits. The single free parameter which is determined in the 
fits is the value of the strong coupling constant. All fit results presented hereafter consider all 
experimental and theoretical uncertainties. The contribution of the uncertainties of the parton 
distributions to the uncertainty of is discussed separately. 

The value of as{Mz) is obtained from a fit to the inclusive jet cross section measured 
double-differentially with the inclusive k_i algorithm. For the result, we use the parton dis- 
tributions from the CTEQ5M1 parameterization [R^ and check the effects of other choices. 



The renormalization scale is chosen to be fir = Et and the factorization scale is set to the fixed 
value of fif = 1/200 GeV (the average Et of the jet sample). The effect on as{Mz) of using a 
different choice for fXj. is studied. 

The studies of the stability of the results include fits to the inclusive jet cross section mea- 
sured with the Aachen jet algorithm, fits to the double-differential dijet cross section 
d^cTdijct/d-ETclQ^ using four different jet algorithms, and fits to other double-differential dijet 
distributions. 



Fits to Single Data Points 

Before carrying out combined fits to groups of data points the consistency of the data is tested 
by performing QCD fits separately to all sixteen single data points of the double-differential 
inclusive jet cross section. 

The fit results are displayed in Fig. [T^ for the four regions of Q^. In each fit a result for 
as{E'r) is extracted which is presented at the average of the corresponding data point. The 
individual results are subsequently evolved to as{Mz)- Combined fits to all four data points in 
the same regions are performed, leading to a combined result of as{Mz) for each region. 
The lower curves in the plots represent the combined fit results and their uncertainties and the 
three upper curves indicate the evolution of the combined result and its uncertainty according to 
the renormalization group equation. The single asiEx) values are consistent with the predicted 
scale dependence of as and all combined as{Mz) results are compatible with each other. The 
results obtained in the different Q"^ regions are (for = Et) 

150 <Q^ < 200 GeV^ 
200 <Q^ < 300 GeV^ 
300 < Q2 < 600 GeV^ 
600 <Q^< 5000 GeV^ : a,(M^) = 0.1188 +°:°°^^ (exp.) t'omll (th.) • (14) 

While the experimental uncertainties are of similar size for all as{Mz) values, the theoretical 
uncertainties shrink slightly towards larger Q^. This is a consequence of the reduced renormal- 
ization scale dependence of the jet cross section at higher Q^. 





{Mz) = 


n 199c: +0.0052 
U.iZZO _o.o054 


(exp. 


\ +0.0060 
1 -0.0062 


(th.) , 




\Mz) = 


n 1 9n9 +0.0044 
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Figure 14: Determination of as from the inclusive jet cross section using the inclusive fc_L 
algorithm at a renormalization scale = . Displayed are the results of the fits to the single 
data points in each Q"^ region at each Ex value (circles) including experimental and theoretical 
uncertainties. The single values are extrapolated to the mass (triangles). A combined fit 
yields a result for as{Mz) (rightmost triangle) for each region. The lower curves represent 
the combined fit results and their uncertainties and the upper curves indicate the prediction of 
the renormalization group equation for their evolution. 

Combined Fit - Central o!.s{Mz) Result 



Having checked that the data are consistent over the whole range of Q"^ and Et combined fits 
are made to groups of data points. To study the E^ dependence of (Et) the four data points of 
the same E^ at different Q"^ are combined and four values of asiEx) are extracted. The results 
are shown in Fig. [T^. The four single values are evolved to as{Mz). A combined fit to all 16 
data points gives x^/n.d.f . = 3.80/15 which is rather small, possibly reflecting a conservative 
estimate of systematic uncertainties. The central result is 



as{Mz) = 0.1186 ± 0.0030 (exp. 



+0.0039 
-0.0045 



(th. 



{fir = E' 



T 



(15) 



in good agreement with the current world average of as{Mz) = 0.1184 ± 0.0031 []45|]. The 
statistical uncertainty of the result is very small (±0.0007). The largest contribution to the 
experimental uncertainty comes from the hadronic energy scale of the LAr calorimeter. The 
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Figure 15: Determination of as from the inclusive jet cross section using the inclusive k±_ algo- 
rithm for the renormalization scale /i^ = Et- The results are shown for each Et value (circles) 
including experimental and theoretical uncertainties. The single values are extrapolated to the 
Z°-mass (triangles). The final result for as{Mz) (rightmost triangle) is obtained in a combined 
fit. The lower curves represent the combined fit result and its uncertainties and the upper curves 
indicate the prediction of the renormalization group equation for its energy evolution. 

theoretical uncertainty includes equal contributions from the uncertainties of the hadronization 
corrections and the renormalization scale dependence. The contribution due to the uncertainty 
of the parton distributions is discussed below. 

Choice of a/Q^ as Renormalization Scale 

Another possible choice of the renormalization scale in the theoretical calculation is the four- 
momentum transfer ^/Q^. Analogous to the procedure applied before, an determination is 
made for the renormalization scale /i^ = ^/Q^■ A combined fit to the 16 data points gives 
X^/n.d.f. = 3.87/15 and a result 

as{Mz) = 0.1227^^034 (exp.) t'ofoll (th.) (/x. = v^) • (16) 

Comparing this result with the one obtained for /i^ = Et in (|l5l), the central value is seen too be 
shifted by +0.0041 and the theoretical uncertainty to have increased substantially. This is due to 
the stronger renormalization scale dependence in the perturbative cross sections for jjr = \fQ^ 
compared to = E^. Within the increased uncertainty contribution from the renormalization 
scale dependence for /i^ = \fQ^ both results are consistent. 
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Figure 16: Dependence of the as{Mz) fit result (for /i^ = Ex) on the parton distributions used 
in the fit. The results are displayed as a function of the as{Mz) value used in the corresponding 
global fit of the parton distributions. The correlation is shown for a comprehensive collection 
of different global fits. 



Using Different Parameterizations of Parton Distributions 

The central fit results are obtained for the parton distributions from the CTEQ5M1 parame- 
terization [01 . The QCD fits are repeated using all parameterizations from recent global 



fits which have been performed in next-to-leading logarithmic accuracy in the MS scheme. 
These include all sets from the fits CTEQ5 [0], Botje99 MRST99 [|§, CTEQ4 gT]], 
MRSR MRSAp [@ and the sets from the gluon uncertainty study by the CTEQ 
coUaborationl^. Many of these fits have provided sets of parton distributions for different as- 
sumptions for as{Mz). Using these sets of parton distributions, the dependence of our results 
on the initially assumed as{Mz) is studied. 

The as{Mz) results obtained for the different parton distributions are shown in Fig. [T^ as a 
function of the as{Mz) value used in the corresponding global fit. The range of the variations 
of the result is small and no significant correlation to the initially assumed as{Mz) is seen. The 
largest deviations from the central result given in (|T5]) are obtained with the MRSR3 parameter- 
ization (+0.0031) and for the set MRST99(gi) (-0.0022). Using the central parameterizations 
from the most recent analyses, results of a, (M^) = 0.1179 for MRST99 and ^^(A/^) = 0.1186 
for Botje99 are obtained which are very close to the result obtained for CTEQ5M1. 

''The parameterization from GRV98 [|T| can not be used since the parameterizations of the charm and the 
bottom quark densities are not provided. 
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Uncertainties in the Parton Distributions 



A determination of the uncertainties of parton density functions (pdfs) has only recently be- 



come available This makes it possible to propagate these uncertainties into the predictions 
of physical quantities. While earlier attempts were restricted to a limited number of variations 
of single parton flavors ^ the fit performed by Botje [ p^ does not only provide para- 



meterizations of the central results, but also the covariance matrix Vij of the 28 fit parameters 
Pi used, including the statistical and experimental systematic uncertainties. In addition further 
systematic studies were performed in by repeating the fit under different physical assump- 
tions; the corresponding deviations are, however, not included in the covariance matrix, but 
presented as single results. The combined information is used to determine the uncertainty of 
the as{Mz) fit result by computing the contributions from the covariance matrix and the single 
systematic studies and add their contributions in quadrature. The uncertainty from the parton 
density functions is then given by 



dpi '•^ dpj 
±0.0019 (pdf: Stat. & exp. 

0.0033 
0.0023 



/ J] ((Aa.(M^))rf 



h0.0027 
-0.0013 



(pdf: fit syst.) 



+0.0033 ^p^f) _ 



(17) 



The largest single contribution comes from the factorization scale dependence which accounts 
for lo;ooo3 Acis{Mz)- In fact, the uncertainty from the parton density functions, determined 
using this procedure is slightly larger than the spread observed in Fig. |16[ The value from ( [T7| ) 
is taken as the uncertainty of our as{Mz) result due to the parton density functions. The final 
result is then 



as{Mz) = 0.1186 ± 0.0030 (exp. 



+0.0039 (i-v. \ +0.0033 f^^f\ 
-0.0045 K^'-'-J -0.0023 IP'^^J 



{fir = Et) ■ 



(18) 



Fits to Other Observables 

To test the stability of the central fit result the same QCD fits are made to some of the other jet 
distributions measured. Included are fits to the differential inclusive jet cross section 
d^cTjct/dE'TdQ^ using the Aachen algorithm and the double-differential dijet cross section as 
a function of various variables for all four jet algorithms mentioned. For the latter the renor- 
malization scale is chosen to be fir = Et- The results of as{Mz) from these fits including 
experimental, theoretical and the pdf uncertainties are displayed in Fig. |T7|. All as{Mz) values 
are in good agreement with each other, with the central fit result given in (|18]) and with the 
current world average value. The results for the exclusive jet algorithms have larger theoretical 
uncertainties due to the larger hadronization corrections. 



5.4 Determination of the Gluon and the Quark Densities in the Proton 

The measurement of described in section |5]^ depends on external knowledge of the parton 
content of the proton, and in particular on the uncertainty in the pdfs of the proton. The deter- 
mined value of Us is found to be consistent with measurements in which no initial state hadrons 
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a from jet production in DIS 
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Figure 17: Comparison of as{Mz) results from fits to different double-differential jet distribu- 
tions. 



are involved, for example in e+e^ annihilation to hadrons [Q. The validity of pQCD at NLO in 
jet production in DIS is thereby demonstrated unequivocally to within the accuracy with which 
the strong coupling constant is known. 

It is therefore appropriate to pursue a determination of the parton density functions of the 
proton in NLO pQCD using measurements of jet production cross sections in DIS, assuming 
the value of the strong coupling constant ag from external measurements. Such a determination 
is important for two reasons. First the measurement is in principle sensitive directly to both 
quark and gluon content in the proton, in contrast with studies of the evolution in xbj and 
of the proton structure function F2 where there is only direct sensitivity to the quark content. 
Second the range in the fractional momentum variable ^ covered by a measurement using jets 
is different from that attained with F2 measurements. 

In the second step of the QCD analysis the sensitivity of the jet cross sections to the gluon 
density in the proton is exploited. The dijet cross section as a function of ^ is directly sensitive 
to the gluon density at a; = ^. The inclusion of the inclusive jet cross section as a function of 
Et maximizes the accessible range in x. Data from a recent measurement of the inclusive DIS 
cross section pE| ] give strong, direct constraints on the quark density xA(x). Furthermore the 
strong coupling constant is set to the world average value as{Mz) = 0.1184 ± 0.0031 [p5|]. 

In the central fit the dijet cross section daj^- Jd^dQ"^ at 150 < < 5 000GeV^ the 
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Figure 18: The gluon density xG{x) in the proton, determined in a combined QCD fit to the 
inclusive DIS cross section, the inclusive jet cross section and the dijet cross section. The jet 
cross sections are measured using the inclusive jet algorithm. The error band includes the 
experimental and the theoretical uncertainties as well as the uncertainty of as{Mz)- 

inclusive jet cross section da'^^^/dETdQ'^ at 150 < < 5 000 GeV^ and the reduced inclusive 
DIS cross section 5-(xBj, Q^) from [|8]] in the range 150 < < 1 000 GeV^ (0.032 < xbj < 
0.65) are used. The gluon density and the x/\{x) quark density are parameterized by 

xP{x) = Ax^ {l-xf il + dx) (19) 

where xP{x) stands for xG{x) or xA(x). 

The gluon density is determined in the range 0.01 < x < 0.1 at the factorization scale 
/i/ = V200 GeV with x^/n.d.f. = 61.16/105. The result is shown in Fig. |1B|. Displayed is 
the error band, including all experimental and theoretical uncertainties and the uncertainty from 
the value of the world average value of as(M^). The result is seen to be in good agreement 
with results from recent global data analyses. The integral of the gluon density over the range 
0.01 < X < 0.1 has been determined to be 

/•O.l 

/ dxxG(x,/iJ = 200GeV2) = 0.229 l[5:|]^J(tot.) , (20) 
io.oi 

= 0.229 l°:°l^(exp.)l»h.)l°:°}^(Aa.). 
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This means that at the scale i^f = V200 GeV gluons with a momentum fraction in the range 
0.01 < X < 0.1 carry 23% of the total proton momentum. This result is in good agreement with 
the results from global fits for which the integral has the values 

CTEQ5M1: 0.227, MRST99: 0.232, GRV98HO: 0.235, Botje99: 0.227. (21) 

The quark density xA{x, fxj = 200 GeV^) determined in this fit is also close to results from 
global fits. To test the stability of the results various cross checks have been performed 

(a) Different parameterizations of the parton densities: 

3 parameters xP{x) = A x'' {I — xY 
5 parameters (I) xP{x) = A x^ {1 — xY {1 + dx^) 
5 parameters (II) xP{x) = A x^ {1 — xY {1 + dy/x + ex) 

While the central result has been obtained using the 4-parameter ansatz in ([T9|), the gluon 
density is unchanged when using other parameterizations and the quark density is stable 
if at least four parameters are used. 

(b) Fits to subsets of the data: The fit has been applied to two subsamples of the data with 
Q"^ < 300 GeV^ and > 300 GeV^ and both, the gluon and the quark resuks are 
unchanged. 

(c) Fits to other jet distributions: The fits have been repeated using other jet distributions 
measured with the inclusive algorithm and also to jet distributions measured with 
other jet algorithms. In all cases the results are consistent with each other. 

5.5 Simultaneous Determination of ck^ and the Proton pdfs 

In the above, as or the gluon density are extracted using external knowledge for the other. A 
more independent test of pQCD can be made in a simultaneous determination of both quantities. 
Such a determination has been performed by fitting the parton densities and as{Mz) using the 
same data sets as in the previous section, the inclusive DIS cross section, the inclusive jet cross 
section (Pa^ct/dErdQ'^ and the dijet cross section d^cTdijct/d^dQ^ (again measured with the in- 
clusive k± algorithm). The gluon and the quark distributions are parameterized according to the 
4-parameter formula in ([T9|). The simultaneous fit yields x^/n.d.f. = 61. 19/104 andaresultfor 
the quark distributions identical to that which is obtained in the fit with a constrained as{Mz). 
The results of this simultaneous fit are displayed in Fig. [1^ as a correlation plot between (M^) 
and the gluon density evaluated at four different values of x = 0.01, 0.02, 0.04, 0.1 which lie in 
the range where the jet cross sections are sensitive. The central fit result is indicated by the full 
marker and the error ellipse is the contour along which the of the fit is by one larger than the 
minimum (including experimental and theoretical uncertainties). The ellipticity of the contours 
indicate that the data included in this analysis are very sensitive to the product ■ xg{x) but do 
not allow a determination of both parameters simultaneously with useful precision. 
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0.1 0.11 0.12 0.13 0.14 0.15 



as(Mz) 



Figure 19: The correlation of the fit results for a^^Mz) and the gluon density at four different 
values of x, determined in a simultaneous QCD fit to the inclusive DIS cross section, the in- 
clusive jet cross section and the dijet cross section. The jet cross sections are measured using 
the inclusive k± jet algorithm. The central fit result is indicated by the full marker. The error 
ellipses include the experimental and the theoretical uncertainties. 



Also included in Fig. [T^ are the results from global fits. All of these results are within 
the error ellipses except for GRV98 (at x < 0.04) which uses a relatively small value of 
asiMz) = 0.114. 

The stability of the results in Fig. [1^ has been tested in a similar way as already described 
in section Fits have been performed excluding either the low (< 200 GeV^) or the high 
(> 600 GeV^) Q'^ data. Although the fits give consistent central results, the high Q'^ data are 
needed to achieve a stable determination of the contour of the error ellipsoid. 



6 Summary 

Jet production has been studied in the Breit frame in deep-inelastic positron-proton collisions 
at a center-of-mass energy of \/s = 300 GeV. In the range of four-momentum transfers 
5 < < 15 000 GeV^ and transverse jet energies 7 < < QO GeV dijet and inclusive 
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jet cross sections have been measured as a function of various variables using k± and angu- 
lar ordered jet clustering algorithms. Perturbative QCD in next-to-leading order in gives 
a good description of all observables for which next-to-leading order corrections are not too 
large and for which hadronization corrections are small. For those observables with moderately 
large hadronization corrections the deviations between data and the perturbative calculations are 
always consistent with the size of the hadronization corrections as predicted by phenomenolog- 
ical models. Only at Q"^ < 10 GeV^ do the theoretical predictions fail to describe the size of 
the measured jet cross sections. In this region, however, the NLO corrections are large, with 
/c-factors above two indicating that NLO calculations are not reliable and that it is likely that the 
perturbative predictions receive large contributions from higher orders in which can account 
for the observed difference. 

QCD analyses of the data have been performed in the region of > 150 GeV^, where 
NLO calculations are reliable, using the inclusive k± jet algorithm for which hadronization 
corrections are smallest. In a first step as has been determined in a fit to the inclusive jet cross 
section as a function of the transverse jet energy. Here the knowledge of the parton density 
functions of the proton is taken from the results of global fits. The observed Et dependence of 
as is consistent with the prediction of the renormalization group equation and a combined fit to 
the data yields 

as{Mz) = 0.1186 ± 0.0059 (tot.) 

= 0.1186 ± 0.0030 (exp.) t'oS (th.) (pdf) . 

This result is seen to be stable when the fit is performed to a variety of jet distributions measured 
with different jet algorithms. 

Including HI data on the inclusive neutral current DIS cross section, the jet data have 
been used for a consistent determination of the gluon density in the proton together with the 
quark densities. Setting as to the world average value within its uncertainty of as{Mz) = 
0.1184 ± 0.0031 the gluon density is determined in the range of momentum fractions 0.01 < 
X < 0.1 at a factorization scale of the order of the transverse jet energies i^f = -\/200 GeV in 
the MS scheme. The integral over the range 0.01 < x < 0.1 is determined to be 

^•0.1 

/ dxx G{x, fj.} = 200 GeV^) = 0.229 taogoCtot.) , 
io.oi 

= 0.229 +&xp.)l°:°^?(th.)l°:°l^(Aa,). 

This result, as well as the differential distribution in x, are in good agreement with the results 
obtained in global fits. 

Finally as and the gluon density in the proton have been determined simultaneously using 
data with direct sensitivity to both. The results and their uncertainties show a large anticorrela- 
tion. Here the single results of ag and the gluon density have relatively large uncertainties, but 
the strong anticorrelation of the combined result clearly demonstrates the high sensitivity of the 
jet data to both. 
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Tables of Experimental Results 



In the following those jet cross sections are listed which have been used in the QCD analyses 
to obtain the central results. The numbers of further distributions can be found in or are 
available on request from the HI collaboration. 



The Inclusive Jet Cross Section d^crjet/(d^^T dQ^) 



bin number 


corresponding range 


1 


150 < < 200 GeV^ 


2 


200 < Q2 < 300 GeV2 


3 


300 < Q2 < 600 GoV^ 


4 


600 < Q2 < 5000 GeV^ 



bin letter 


corresponding Et range 


a 
b 



d 


7 < Et < llGcV 
U< Et < 18GeV 
18 < Et < 30 GeV 
30 < Et < 50 GeV 



the inclusive jet cross section 



dVjet/(dETdQ2) 



— inclusive k± jet algorithm 





















single contributions to correlated uncertainty 






bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model dep. 


positron 


positron 


LAr hadr. 


hadroniz. 


No. 


section 


uncert. 


uncertainty 


uncertainty 


uncertainty 


detector corr. 


energy scale 


polar angle 


energy scale 


correct. 




(in pb) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(percent) 
























i a 




±S.S 


9.6 


-9.5 


7.6 


-7.4 


5.9 


-5.9 


±4.7 


O.S 


1.7 


2.6 


-2.6 


-7.g 


1 b 


26.084 


± 4.4 


16.3 


-13.9 


13.6 


-11.3 


9.0 


-8.0 


± 6.3 


0.8 


1.6 


6.0 


-4.4 


-5.0 


1 c 


5.819 


± 9.2 


13.8 


-15.6 


12.3 


-13.8 


6.3 


-7.2 


± 4.3 


1.0 


1.7 


3.9 


-5.3 


-4.5 


1 d 


0.719 


± 27.8 


34.0 


-35.4 


32.4 


-33.4 


10.4 


-11.6 


± 3.4 


1.9 


2.7 


9.2 


-10.5 


-4.7 


2a 


62.256 


± 2.6 


7.6 


-7.2 


6.4 


-6.1 


4.1 


-3.9 


± 2.3 


0.2 


1.0 


2.8 


-2.5 


-8.1 


2b 


29.802 


± 3.7 


12.2 


-12.7 


10.3 


-10.8 


6.5 


-6.7 


±4.0 


0.7 


0.2 


4.9 


-5.2 


-4.5 


2c 


6.989 


± 7.6 


16.8 


-14.8 


14.3 


-12.4 


8.9 


-8.1 


± 5.7 


0.3 


1.9 


5.3 


-3.8 


-4.6 


2d 


0.994 


± 20.0 


28.7 


-33.5 


26.0 


-29.9 


12.1 


-15.0 


± 8.6 


2.7 


1.9 


7.7 


-11.8 


-4.9 


3a 


61.577 


± 2.7 


5.8 


-6.1 


4.9 


-5.3 


3.0 


-3.1 


± 1.2 


0.7 


0.8 


2.0 


-2.2 


-8.1 


3b 


35.010 


± 3.5 


11.9 


-9.9 


10.1 


-8.2 


6.3 


-5.5 


± 3.5 


0.5 


1.4 


4.8 


-3.6 


-4.0 


3c 


9.644 


± 6.6 


15.7 


-17.9 


13.3 


-15.3 


8.4 


-9.4 


± 4.2 


1.6 


3.7 


5.9 


-7.2 


-4.5 


3d 


1.362 


± 20.0 


36.3 


-29.7 


32.1 


-26.4 


17.1 


-13.7 


±11.4 


1.1 


1.6 


12.5 


-7.1 


-4.6 


4a 


46.515 


± 3.1 


7.4 


-6.8 


6.3 


-5.8 


3.8 


-3.5 


± 0.7 


0.0 


1.4 


3.1 


-2.7 


-9.0 


4b 


26.409 


±4.1 


8.8 


-8.8 


7.6 


-7.6 


4.4 


-4.4 


± 1.6 


0.6 


1.5 


3.5 


-3.5 


-4.0 


4c 


11.288 


± 6.0 


11.5 


-11.6 


10.3 


-10.3 


5.3 


-5.4 


±2.4 


0.3 


0.5 


4.4 


-4.6 


-3.2 


4d 


1.993 


± 15.1 


27.2 


-23.0 


24.6 


-21.2 


11.5 


-8.9 


± 1.2 


0.9 


1.9 


11.1 


-8.4 


-3.2 



Table 2: Results of the inclusive jet cross section measurement using the inclusive k± jet al- 



gorithm. The listing includes all experimental uncertainties (as described in section |J!4| ) which 
are here separated into the correlated and the uncorrelated part. Since the interpretation of the 
results (as e.g. in a QCD analysis) does not require the knowledge of the separate contributions 
to the uncorrelated part of the uncertainties, only the total uncorrelated uncertainty is presented 
while the single contributions to the correlated uncertainty are listed in extra columns for all 
sources. The uncertainty from the hadronic energy scale of the Liquid Argon calorimeter is 
quoted asymmetrically. The left (right) value corresponds to an increase (decrease) of the cal- 
ibration constants. The uncertainties of the positron energy and the positron polar angle are 
defined to be symmetric by taking the maximum of the upwards and downwards deviations. 
The signs are quoted for a positive variation of the corresponding source. Note that only the 
correlated contribution from these sources is listed. As described in section WM some of these 



sources contribute also to the uncorrelated uncertainty. The latter contribution is already con- 
tained in the (quadratic) sum of all uncorrelated uncertainties. The total correlated uncertainty 
includes also the contribution of ±1.5% from the uncertainty in the determination of the lumi- 
nosity. In the right column we have also included the size of the hadronization corrections as 
determined by the procedure described in section 
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The Dijet Cross Section d?(TAi3et/{di dQ"^) 



bin No. 


corresponding range 


5 range 


1 a 


5 < < lOGeV^ 


0.004 < 5 < 0.01 


1 b 




0.01 < 5 < 0.025 


1 c 




0.025 < 5 < 0.05 


1 d 




0.05 < 5 < 0.1 


2a 


10 < < 20 GeV^ 


0.004 < 5 < 0.01 


2b 




0.01 < 5 < 0.025 


2c 




0.025 < 5 < 0.05 


2d 




0.05 < 5 < 0.1 


3a 


20 < < 35 GeV^ 


0.004 < 5 < 0.01 


3b 




0.01 <i< 0.025 


3c 




0.025 < ? < 0.05 


3d 




0.05 < 5 < 0.1 


4a 


35 < < 70 GeV^ 


0.004 < 5 < 0.01 


4b 




0.01 <i< 0.025 


4c 




0.025 < ^ < 0.05 


4d 




0.05 < ? < 0.1 



bin No. 


corresponding range 


5 range 


5a 


150 < < 200 GeV^ 


0.009 < C < 0.017 


5b 




0.017 <^< 0.025 


5c 




0.025 < 5 < 0.035 


5d 




0.035 <i< 0.05 


5e 




0.05 < ^ < 0.12 


6a 


200 < < 300 GeV^ 


0.01 <£,< 0.02 


6b 




0.02 <^< 0.03 


6c 




0.03 < 5 < 0.04 


6d 




0.04 < 5 < 0.06 


6e 




0.06 < C < 0.15 


7a 


300 < < 600 GeV^ 


0.015 < C < 0.025 


7b 




0.025 < 5 < 0.035 


7c 




0.035 < C < 0.045 


7d 




0.045 < ^ < 0.065 


7e 




0.065 <i< 0.18 


8 a 


600 <Q'^ < 5000 GeV^ 


0.025 < ^ < 0.045 


8b 




0.045 <i< 0.065 


8c 




0.065 < 5 < 0.1 


8d 




0.1 < 0.3 





the dijet cross section 


d'^o-di^ct/id^dQ'^) — inclusive A; 


j_ jet algorithm 






















single contributions to correlated uncertainty 






bin 


cross 


statistical 


total 


uncorrelated 


correlated 


model dep. 


positron 


positron 


LAr hadr. 


hadroniz. 


No. 


section 


uncert. 


uncertainty 


uncertainty 


uncertainty 


detector corr. 


energy scale 


polar angle 


energy scale 


correct. 




(in pb) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(in percent) 


(percent) 


1 a 


19.34 


± 6.8 


10.6 


-10.0 


8.5 


-8.5 


6.3 


-5.3 


± 2.4 


-0.2 


-3.6 


4.9 


-3.0 


-6.1 


1 b 


83.84 


± 3.8 


10.8 


-9.2 


8.4 


-6.3 


6.8 


-6.7 


± 1.1 


2.7 


4.8 


3.7 


-3.6 


-9.3 


1 c 


47.96 


± 5.2 


6.3 


-10.1 


4.8 


-8.2 


4.2 


-5.9 


± 3.8 


0.6 


2.7 


2.7 


-3.5 


-5.9 


1 d 


18.72 


± 7.4 


16.8 


-8.6 


13.6 


-6.9 


9.9 


-5.2 


± 6.2 


2.6 


2.1 


6.8 


-4.0 


-9.3 


2a 


14.67 


± 7.1 


8.6 


-9.2 


6.9 


-7.6 


5.1 


-5.2 


± 2.2 


1.6 


-2.3 


4.0 


-3.8 


-13.2 


2b 


66.71 


± 3.6 


8.0 


-9.9 


6.6 


-8.7 


4.5 


-4.7 


± 2.0 


1.0 


-1.6 


3.8 


-3.8 


-6.4 


2c 


39.42 


± 5.0 


12.2 


-7.8 


10.2 


-5.6 


6.6 


-5.4 


± 3.2 


1.1 


3.2 


5.4 


-2.7 


-5.6 


2d 


14.58 


± 7.8 


11.2 


-13.3 


8.8 


-10.8 


6.9 


-7.7 


± 3.7 


-1.7 


4.0 


5.1 


-5.2 


-11.4 


3a 


9.45 


± 8.8 


15.7 


-11.9 


13.4 


-10.6 


8.2 


-5.6 


± 1.7 


2.4 


1.1 


7.6 


-4.6 


-12.7 


3b 


49.42 


± 4.0 


7.9 


-9.2 


6.5 


-8.1 


4.5 


-4.4 


± 1.1 


0.8 


1.7 


3.8 


-3.8 


-7.0 


3c 


28.83 


± 5.3 


7.1 


-11.4 


5.4 


-9.9 


4.6 


-5.7 


± 0.3 


-1.8 


-2.4 


3.1 


-4.8 


-6.6 


3d 


10.90 


± 9.6 


12.4 


-12.4 


10.2 


-9.8 


7.1 


-7.6 


± 4.7 


3.1 


1.8 


5.9 


-4.8 


4.2 


4a 


6.46 


± 10.6 


16.7 


-9.3 


15.0 


-7.6 


7.4 


-5.3 


±1.5 


1.9 


1.9 


6.7 


-4.3 


-3.0 


4b 


47.92 


± 3.8 


9.5 


-7.0 


8.2 


-5.7 


4.9 


-4.1 


± 1.0 


-1.3 


2.1 


3.8 


-3.3 


-4.7 


4c 


31.09 


±4.9 


6.7 


-9.8 


5.4 


-8.3 


3.8 


-5.3 


± 2.4 


-1.1 


1.1 


3.1 


-4.4 


-5.0 


4d 


12.06 


± 8.0 


13.8 


-9.0 


11.9 


-7.5 


7.1 


-5.0 


± 0.4 


0.9 


2.4 


6.5 


-4.3 


-11.5 


5a 


4.147 


±11.0 


24.7 


-27.7 


19.9 


-22.9 


14.7 


-15.6 


± 13.7 


1.9 


1.8 


4.3 


-6.9 


-6.0 


5b 


6.272 


± 9.1 


14.9 


-13.8 


13.2 


-12.2 


6.9 


-6.4 


± 3.4 


2.1 


2.8 


4.7 


-3.8 


-5.3 


5c 


6.544 


± 9.8 


12.3 


-13.2 


11.2 


-12.0 


5.1 


-5.6 


± 3.1 


3.0 


1.2 


1.8 


-3.0 


-5.7 


5d 


5.059 


± 10.9 


14.3 


-15.3 


13.2 


-14.1 


5.3 


-5.9 


± 1.9 


2.3 


2.2 


3.5 


-4.4 


-5.4 


5e 


4.800 


±11.3 


18.4 


-17.8 


15.1 


-14.5 


10.6 


-10.3 


± 7.9 


3.3 


5.3 


2.9 


-1.7 


-6.2 


6a 


6.324 


± 8.8 


13.7 


-18.1 


12.0 


-16.0 


6.5 


-8.6 


± 3.8 


3.6 


1.2 


3.4 


-6.5 


-6.1 


6b 


7.309 


± 7.8 


10.9 


-12.5 


10.0 


-11.3 


4.4 


-5.3 


± 2.6 


0.3 


1.0 


3.1 


-4.3 


-4.8 


6c 


6.023 


± 8.6 


11.7 


-11.8 


10.5 


-10.5 


5.2 


-5.5 


± 1.1 


1.5 


3.7 


2.8 


-3.2 


-5.0 


6d 


5.512 


± 8.9 


14.9 


-11.9 


13.4 


-11.0 


6.6 


-4.6 


± 2.6 


1.6 


1.9 


5.3 


-2.4 


-4.8 


6e 


4.186 


± 10.3 


17.8 


-19.3 


14.7 


-16.1 


10.1 


-10.8 


± 8.3 


3.1 


3.6 


3.0 


-4.8 


-7.5 


7a 


5.997 


± 9.3 


12.7 


-15.4 


11.9 


-14.4 


4.5 


-5.6 


± 1.5 


0.9 


1.0 


3.8 


-5.1 


-5.4 


7b 


7.006 


± 8.3 


11.3 


-12.4 


10.1 


-11.1 


5.0 


-5.6 


±3.5 


0.4 


2.1 


2.4 


-3.4 


-4.8 


7c 


6.104 


± 9.0 


16.0 


-11.9 


14.5 


-11.1 


6.9 


-4.5 


± 1.2 


2.2 


1.0 


6.1 


-3.3 


-4.5 


7d 


7.249 


± 8.1 


10.1 


-10.4 


9.0 


-9.2 


4.7 


-4.9 


± 1.9 


0.6 


3.8 


1.2 


-1.9 


-6.2 


7e 


6.082 


± 8.8 


14.7 


-14.3 


13.1 


-12.7 


6.7 


-6.5 


± 3.0 


1.2 


2.4 


5.1 


-4.9 


-5.8 


8a 


6.077 


± 9.4 


15.0 


-12.6 


13.6 


-11.6 


6.3 


-5.0 


± 2.3 


0.6 


2.0 


5.2 


-3.5 


-6.2 


8b 


6.759 


± 8.5 


13.3 


-11.9 


11.7 


-10.5 


6.3 


-5.6 


± 3.8 


2.3 


1.9 


3.8 


-2.4 


-5.6 


8c 


8.305 


± 8.0 


11.0 


-10.2 


10.2 


-9.5 


4.1 


-3.5 


± 0.6 


0.9 


1.0 


3.6 


-2.8 


-4.8 


8d 


7.520 


± 8.2 


10.2 


-11.2 


9.6 


-10.4 


3.3 


-4.1 


± 0.4 


1.3 


0.3 


2.7 


-3.5 


-5.7 



Table 3: Results of the dijet cross section measurement at high Q'^ using the inclusive k± jet 
algorithm. The presentation is as in table ^ 
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Results of the QCD Analysis 

Et Dependence of ck^ 



Et dependence 0^ as {Et) (fir ^ Et) 
inclusive jet cross section — inclusive k± algorithm 


average Et of data point: 


V70 GcV 


V200 GcV 


V500 GcV 


VI 500 GcV 


asiET) = 


0.1940 


0.1636 


0.1579 


0.1507 


total uncertainty 


+U.UibV 
-0.0145 


+0.0145 
-0.0136 


+0.0162 
-0.0167 


+0.0308 
-0.0315 


exp. 

theor. 

pdf 


+0.0082 
-0.0081 
+0.0097 
-0.0103 
+0.0092 
-0.0061 


+0.0106 
-0.0104 
+0.0077 
-0.0081 
+0.0040 
-0.0033 


+0.0124 
-0.0123 
+0.0091 
-0.0086 
+0.0051 
-0.0072 


+0.0264 
-0.0244 
+0.0113 
-0.0088 
+0.0110 
-0.0179 


as{Mz) - 


0.1211 


0.1174 


0.1227 


0.1292 


total uncertainty 


+0.0058 
-0.0056 


+0.0068 
-0.0070 


+0.0094 
-0.0101 


+0.0218 
-0.0233 


exp. 

theor. 

pdf 


+U.0031 
-0.0031 
+0.0035 
-0.0040 
+0.0035 
-0.0023 


+0.0052 
-0.0053 
+0.0038 
-0.0043 
+0.0020 
-0.0017 


+0.0U72 
-0.0074 
+0.0053 
-0.0054 
+0.0030 
-0.0042 


+0.0187 
-0.0182 
+0.0079 
-0.0067 
+0.0078 
-0.0129 



Table 4: The ag results from the fits presented in section Displayed are the fit results of 
asiEr) at different Et (top) and the corresponding values extrapolated to fir = Mz (bottom) 
together with the different contributions to the uncertainty. 



The Gluon Density in the Proton 



The Gluon Density in the Proton at a</ = \/200GcV 
determined for as{Mz) = 0.1184 ± 0.0031 
parameterized by xG{x) = Ax'' (1 - x)" (1 + dx) in 0.01 < a; < 0.1 
central result: A=0.503 ; b = -0.5935 ; c = 4.70 ; d = -0.55 



io(a;) 


xG{x) ^ 


-2.0 


7.35 


-1.9 


6 32 +1"'' 

"•'^'^ -0.98 


-1.8 


r 42 +0-82 
"^•^^ -0.77 


-1.7 


4 ficj +0.65 
^•O-J -0.63 


-1.6 


o OQ +0.54 
•^■^•^ -0.52 


-1.5 


q qi +0.45 
o.Oi _o.43 


-1.4 


2 76 

^- -0.36 


-1.3 


9 97 +0.31 
-0.30 


-1.2 


1 04 +0.26 


-1.1 


1 47 +0-22 
-0.20 


-1.0 


1 14+019 
^•-L^ -0.17 



exp. 

+0.93 
-0.94 

+0.64 
-0.65 
+0.45 
-0.45 
+0.32 
-0.32 
+0.26 
-0.25 
+0.23 
-0.21 
+0.21 
-0.20 
+0.19 
-0.18 
+0.17 
-0.16 
+0.15 
-0.14 
+0.13 
-0.13 



theor. from Aas{Mz) 
+0.57 



+0.77 

-0.71 
+0.65 
-0.63 
+0.51 
-0.54 
+0.44 
-0.46 
+0.36 
-0.38 
+0.28 
-0.31 
+0.22 
-0.25 
+0.18 
-0.19 
+0.13 
-0.14 
+0.12 
-0.11 
+0.08 
-0.10 



-0.43 
+0.49 
-0.37 
+0.42 
-0.32 
+0.36 
-0.28 
+0.30 
-0.24 
+0.25 
-0.20 
+0.21 
-0.17 
+0.17 
-0.14 
+0.14 
-0.12 
+0.11 
-0.09 
+0.08 
-0.07 



Table 5: The gluon density in the proton from the fit in section |5]^. Displayed are the central 
results and the total uncertainties of the gluon density at eleven values of x in the interval 
0.01 < X < 0.1 together with the different contributions to the uncertainty. Also shown are the 
parameters A, b, c, d of the central result. 
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